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We study the free energy of the pure glue QCD string with a torus target space and the 
gauge groups SU(N) and (chiral) U(N). It is highly constrained by a strong/weak gauge 
coupling duality which results in modular covariance. The string free energy is computed 
exactly in terms of modular forms for worldsheet genera 1-8. It has a surprisingly mild 
singularity in the weak gauge coupling/small area limit. 
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1 Introduction 

The idea of a string theoretic formulation of QCD is as tantalizing today as it was twenty 
years ago. Despite its age and elusiveness, the promise of a description of the phenomena 
of strongly coupled gauge theory in term of strings is compelling. The many problems 
encountered in trying to implement this idea have shown that the formulation must differ 
substantially from the critical string. It must not include spacetime gravity, for instance. As 
a result, there are no promising models at this time. 

The hope for a QCD string has been rejuvenated recently by significant progress un- 
derstanding QCD in two dimensions. Pure glue QCD 2 has been solved exactly |l|, and a 
great deal of work has gone into constructing a string theory from its large N expansion 
|2|. ^ |J . The string action that has resulted is a perturbed topological sigma model coupled 
to topological gravity || ||, 0, ||. Perhaps a simpler formulation of the topological string is 
possible, but we are in a novel situation of having a string theory that has been shown to be 
equivalent to QCD, albeit in two dimensions. 

The ultimate goal is to use the string theory of QCD 2 to construct, directly or indirectly, 
a QCD 4 string action. The direct extension of the sigma model to four dimensions would 
presumably yield some version of topological Yang- Mills, although this is a non-trivial, open 
problem. A further perturbation would be necessary to get dynamical (pure glue) QCD. 

If the direct approach proves to be intractable, we may still make progress toward a 
QCD 4 string by identifying properties of the string that are independent of dimension. To 
that end, we will examine the QCD string partition function in detail, particularly focusing 
on QCD on the torus which most resembles critical string theory. The large N expansion of 
the QCD free energy has been given in terms of two different group theoretic sums |l|, |J, the 
free energy for free fermions and the equivalent Jevicki-Sakita bosons || |TU| and, of course, 
the free energy of the topological string itself ||. The various formulations readily reveal 
different aspects of the free energy, but none is completely explicit. We will show how the 
structure of the Jevicki-Sakita expressions and the simple algebra of the heat kernel sum 
combine to let us compute the free energy efficiently. The free energy will be calculated 
exactly in terms of modular forms up to genus 8. 

It is very rare to have an exact expression to the eighth order of string perturbation 
theory. In those cases where such an expression has been found, as in matrix models, it is 
possible to continue to all orders. In fact, there is a differential equation relating the free 
energy at a given order to that at lower orders. If such an equation were known for the QCD 
string, then there would be no need to display the horrendous expression for genus 8. But 
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the obvious candidates (slight generalizations of the holomorphic anomaly equation |J) fail. 
The absence of boundary contributions to the (chiral) U(N) free energy suggests an even 
simpler structure-that a handle creation operator exists. We find something akin to one, 
but it couples to the infinitely many deformations of QCD. A simple equation in terms of 
the Kahler modulus alone would be much more powerful. It is plausible, but it has not been 
found. So we will list the exact expressions up to genus 8 and point out some surprising 
features that emerge. 

2 QCD2 and Its String Expansion 

This section presents an overview of some of the salient aspects of QCD2. We give a brief 
review of the heat kernel partition function and Gross's large N expansion of it. Next we 
discuss the relationship between SU(N) and U(N) 2D Yang-Mills theory, and we show that 
the SU(N) partition function is easily computed from that of chiral U(N). 

Pure glue Yang-Mills theory is exactly solvable in two dimensions. This is largely due to 
the absence of propagating gluons, as only global degrees of freedom survive gauge fixing. 
There are no transverse gluons to propagate. The partition function is 



Z = J[DA,}e~£^ F ^ d2x (2.1) 



which may be calculated on any Riemann surface E. It only depends on the scaled gauge 
coupling A (= QqcdN), the area A of E, the topology (genus G) of E and the gauge group. 
We will consider Riemann surfaces with no boundaries. The gauge group will be either 
SU(N) or U(N), with a large N in order to get the string expansion in g st = 1/N. 

There is a remarkable solution of QCD2 due to Migdal and Rusakov [|]] . The heat kernel 
lattice action reproduces fl2.1|) in the continuum limit, and it has the powerful feature that it 
is renormalization group invariant. This permits a quick solution for the partition function 

Z = £(dim Rf-^ e -^c 2 (R) (22) 

R 

where C^-R) is the second Casimir and dimi? is the dimension of the representation R. 
The sum over irreducible representations may be expressed as a sum over different weights 
(i.e. Young tableaux). A Young tableau has n^ boxes in the k th row and rows decreasing 
in length, r\,\ > n-i > • • • > n^. For SU(N) each n^ is a non- negative integer, whereas for 
U(N) the weights n& may be any integer. The additional irreducible representations are due 



to the 17(1) in U(N) = (SU(N) x U{1))/Z N . The Casimir is 



C 2 {R) = 



N 



rr 



J2n k (n k + N + l-2k)-— for SU(N) 
fc=i 



iV 



(2.3) 



J2Mn k + N + l-2k) 
K. fc=i 



for U(N) 



where n = J2 n k and a particular choice is made for the U(l) charge in order to simplify the 
U(N) Casimir. (In general, there is an extra term ^-, where a is determined by the U(l) 



charge y X' in U(N): a = 1 — NJX'/X.) The dimension of R is given by 



dimi? 



It is a polynomial in N of degree n, 



n 

l<i<j<7V 



rii - rij - i + j 
3 -i 



(2.4) 



dimi? = - J2 N K °Xr(.v) 

crGbn 
d ( °° 



%****, {u-V%*^ 



.k=l 



(2.5) 



The first expression comes from the Frobenius formula for characters X-r(°~) °f the represen- 
tation R of the symmetric group S n [|J. K a is the number of cycles in the permutation a . 
The second expression may be considered the large N expansion of the standard "factors 
over hooks" rule for the dimension (cf. ||11| ). The invariants C^(R) are given by 



N rij 

<?«(*) = EEC*-* 



fc-1 



(2.6) 



,U(N), 



C( fc) (-R) is part of the k th Casimir, C k (R) = k\C (k) (R) + ■■■. For example, C^ '{R) 
nN + 2(7(2) (R) an d n — C(i}(R). The resulting expression for the partition function is 



n.j-rij-i+j 



2-2G \A 



J2 k [n k (n k +N+l-2k)-§n 2 



N-l 



z = e n 

rii>ni+i l<i<j<N 

[fe»>hi+i l<i<j<N J Z=l 

where h k = n k + ~(iV + 1) - fe (cf. 0) and n = J2 n k = Y.h k . 



2G-2 



(2.7) 



The U(N) case is simpler both because of the form of the Casimir and because the sums 
run over all integers nk, including the negative ones. This lets us express the G = 1 partition 
function in terms of elliptic functions, 



A^! 



. e ^f(iv 2 -i) 



7?(t) JV -(^)^(t) JV - 2 ^(2*) 



where £ = i\A/2nN. The corrections come from terms in $(t) N with /ij = /i^ for some % ^ j. 
The Jacobi theta function d is $2 when A" is even and $3 when A" is odd. 

The SU(N) partition function is more complicated. Consider first the U(N) partition 
function sum restricted to SU(N) Young tableaux. It may be expressed in terms of the 
functions 

00 1 1 (AT-3)/2 

V)(*)= E **** = -<?(*) + - E e^ a 

h=(3-N)/2 A Z h=(3-N)/2 / 2 g x 

00 

h=(N-l)/2 

Then the intermediate partition function is 



itfe 2 



^ = _L_^ 3 -3^-3) | V)(i) ^-i - ( V ) V)^ 3 V)(2t) 



J- ^4-fiV3_SA7 2 ^KiV-S f n /.x/V-1 f N— 1 

(TV^l)! 

(2.10) 

There is a related function, Z + , in which the sum is restricted to Young tableaux with fewer 

than N/2 boxes in any column. 



Z+ = ^ e ^(^-D [ut) N/2 - 2 ( ^f ) tf 2 (^- 4 )/ 2 ^ 2 (2t) - ■ ■ ■} (2.11) 

(for A" even) . The large A" expansion of this seemingly ad hoc function will be the focus of 
much of what follows. Finally, to get the SU(N) partition function, the extra term in the 
Casimir must be included, but that is difficult. In any case, these expressions are curious, 



but they are not much help. They work well for small N: 






-7(3=1 

Zj U{2) 



yG=l 

Z t/(3) 
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Zj U(2)/SU 
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Zj U(3)/SU 



7 G=\ 

Zj SU(2) 



§e¥{0 2 (t) 2 -0 2 (2t)} 

|e^ {& 3 (t) 3 - M 3 (t)M2t) + 2^ 3 (3t)} 

i e ¥^3(t)2 + 2^(t)-2^ 3 (2t)-l} 



(2.12) 



\e^ {urn) - 1} 



Unfortunately, ( |2.8|) and ( 2.11 ) are not conducive to large A" expansions, since it is difficult 
to determine if one term dominates the sum. Also, note that even the simpler U(N) partition 
function expressed in terms of elliptic functions is not a modular form, since each term has 
a different weight. We will see below that at each worldsheet genus the G = 1 string free 
energy is almost a modular form, but with a different modulus r = yt. On the other hand, 
the modular weights of the theta functions do determine the small area behavior of the 
partition function. Since ^(t) ~ £~ 5 as t — > 0, the leading term dominates the small area 
limit at finite N. If A" is then taken to infinity, the partition function develops an essential 
singularity at t — 0. It is not clear from (|2.8| ) if one term dominates in the large A^ limit, 
which should be taken first. We will see below that the string partition function has an 
essential singularity at XA = which is the phase transition that occurs at finite coupling 
on the sphere |T2[ . 

The string expansion for Z is a large A^ expansion. This is explained in detail in and 
||, but a brief discussion of the structure will help motivate the ensuing analysis. Gross 
and Taylor have shown how SU(N) representations with relatively small Casimirs (of order 
N) and small dimensions give the leading contribution to the partition function at large N, 
yielding a series that has many properties of closed string perturbation theory. A Young 
tableau with a small number of boxes has a relatively small Casimir, and it makes a leading 
contribution. But for SU(N), a representation R and its complex conjugate R have the 
same Casimir and the same dimension, so they make the same contribution to Z. This 
leads to a natural factorization of representations. The "crural" representations are those 
with no more than A^/2 boxes in any column, and the "anti-chiral" representations are those 



whose complex conjugate is chiral with no column of N/2 boxes.[] (Recall that if R has Cj 
boxes in its j th column, then R has N — Cj boxes in its j th column from the right.) Then 
any representation is expressed uniquely as the Young product of an anti-chiral and a chiral 
representation; i.e. any tableau is a chiral tableau joined to an anti-chiral tableau. 

The physical partition function may be obtained from the chiral U(N) partition function, 
Z + , in which the sum is restricted to chiral SU(N) representations, 



\A, 



Z + = J2(dimR) 2 - 2G e-^ 2C ^e- nA (2.13) 

R+ 

This is a well-defined function whose asymptotic expansion is the holomorphic topological 
string perturbation theory. A = |AA, but it is kept formally independent of A so that the 
extra piece of the SU(N) Casimir may be obtained by differentiation with respect to A. For 
example, the chiral SU(N) partition function is given by 



J SU(N) 

The free energy F = log Z is 



Z+=e^n Z +{A) 



A=±\A 



(2.14) 
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2N 2 

XA \2 

2N 2 > 



+ ■■■ 

(2.15) 
where the Schur polynomials P n are generated by e XkZ = yj P n (xi, ■ ■ • , x n )z n , the polynomi- 
als P n are generated by log(l + XkZ k ) — Y^Pni^i, • ■ ■ , x n )z n , and the primes denote d^. The 
full SU(N) partition function on the torus is 



^ G = 1 )- e (^)^ 1 -^) 2 z+(i 1 )z + (i 2 ) 



. _ 1 (2-16) 



^SU(N) 

using the formula for the Casimir of the Young product of a chiral and an anti-chiral repre- 



1 This is a slight generalization of the definition given by Gross and Taylor M. 



sentation, C 2 (SR) = C 2 (R) + C 2 (S) + 2nn/N. The G = 1 free energies are related by 

XA _ 



F = 2F++Y: Pm fa = iP 2i (0, 2F+", 0, 2F+W, • ■ ■ , 2F + ^\ 
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F +W + 6 (F + ") 2 ] + ^^- [F +(6) + 30 F + " F + ^ + 60(F + ") 3 

F +(8) + 56F+"F + ( 6 ) + 70(F + ^) 2 + 840 ((F+") 4 + (F + ") 2 F + ( 4 )))] + 

(2.17) 
Note that only even derivatives of F + enter ( |2.17| ). As a result F + and F (but not Fg V , N s 
or F(/(tv)) have a simple modular structure independent of the zero point energy. Once the 
chiral U(N) free energy is known, it is straight-forward to calculate Fgy,^ and Fsu(n)- Of 
course, it is trivial to get Fjj, N <. from the full SU(N) free energy, F S u(n), because of its 
modular structure. 

The form of ( [2.16 ) suggests that for G = 1 the U(N) chiral partition function is on the 



same footing as the more physical non-chiral SU(N) partition function. This is remarkable 
since a topological string theory reproducing simple Hurwitz space (i.e. the moduli space of 
chiral U(N), cf. 0) is relatively easy to construct. There are no contributions from the 
boundary of moduli space; the contact terms vanish. A very simple, explicit perturbation 
gives the full theory. This is in striking contrast to the complications of G ^ 1. 

The formula for the dimension of a composite representation is not especially simple, so 
( p.l6j ) and (|2.17|) require unwieldy corrections for G^l. In fact, of all the group invariants, 
only the quadratic Casimir has such a simple decomposition. The higher Casimirs, Ck{SR), 
decompose into a sum of products of the lower Casimirs Ci(R) and Ci(S) with / < k. The 
full partition function of QCD perturbed by C&, may be expressed as a deformation of two 
copies of Z + depending on the couplings of all the lower Cj's. The higher the Casimir, the 
more couplings that must be differentiated. The dimension of R may be expressed in terms 
of the Cj's as well, but in the large N limit it takes infinitely many, so the analog of ( |2.16| ) 
for G^l requires derivatives with respect to infinitely many couplings. 

The perturbed chiral U(N) partition function may be written 

Z+(A,A 2 , ■■■) = ]T(dimi?) 2 - 2G exp ^Ti^i?)/^- 1 ) e' nA . (2.18) 

R+ V k I 

This partition function results from the chiral reduction of a renormalization group invariant 
heat kernel lattice action as before, but it leads to perturbations of the Yang-Mills action 
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F 2 by higher powers of the field strength (F 2 ) k in the continuum. Its large N expansion is 
string-like, since only even powers of 1/N arise. This results from a cancellation in the sum 



over representations. It follows from the definition (|2.6D that Cqa{R) = (— 1) 1 Cn l; -\(S), 



if the columns of .R's Young tableau are the rows of S"s. Since the Casimir Ck(R) is a 
homogeneous polynomial in C \i){R) / N 1 " 1 , the relative minus sign cancels any occurrence of 
an odd power of 1/N arising from the perturbations or the dimension ( |2.5|) . This argument 
immediately extends to the SU(N) case and the full partition function. Of course, there are 
other properties a string perturbation expansion should possess, but these will be verified 



elsewhere 13 



The large N chiral partition function for U(N) Yang-Mills theory on the torus has been 
shown to be the following sum over classes of n-sheeted holomorphic covering maps v n Q 

oo oo 

F+ = E^ 2 " 29 E E ^(lAAl'e^ (2.19) 

g =G n=l [i/ n ]:£ 9 -»E G= i 

where £ 9 (S G ) is the genus g (G) worldsheet (space-time) Riemann surface, and [u n ] is a class 
of n-sheeted covers with i square root branch points. f9„ is the number of holomorphic 
maps divided by a symmetry factor characterizing the moduli space of maps. The maps 
satisfy the Riemann-Hurwitz relation 

2(0-1) = 2n(G-l)+z. (2.20) 



Since i = 2g — 2 for maps to the torus, the factors multiplying the exponential in ( j2.19| ) 



combine into (XA/2N) 29 2 . The QCD partition function on other (G 7^ 1) Riemann surfaces 
is more complicated than (|2. 19 ). Even the SU(N) partition function on the torus is more 



complicated, with contributions coming from the boundary of moduli space where worldsheet 
handles are collapsed to a point in space-time 0. For G 7^ 1 at XA = 0, the free energy is a 
sum of the orbifold Euler characteristics of the moduli spaces 0] . 

We will study QCD on the torus, so there is no contribution to the free energy coming 
from genus worldsheets. The lowest genus worldsheet to contribute is an unbranched single 
cover of the target space, so g m i n — G according to the Riemann-Hurwitz formula. Since 
Fq = the relationship between the free energy and the partition function is simple: 

Z = e F ^exp\Y^(l/N) 29 - 2 F g 

^= 2 x (2.21) 



e Fl 



(l/iV) 2 F 2 + (1/iV) 4 (F 3 + \F 2 ) + • 



We will usually discuss the free energy below, but the partition function is very similar in 
form. For instance, when F + has modular weight 6(7 — 6, so does Z+ apart from the factor 
e F i = fj. This eta function factor gives the partition function an essential singularity at zero 
area/coupling which can cause difficulties (like order of limits problems) that do not arise 
with the free energy. 

3 QCD On The Torus 

In this section we will proceed to study some properties of the G — 1 free energy in detail. 
The free energy on the torus has an interesting structure reminiscent of critical string theory, 
but not shared with QCD on other Riemann surfaces. It transforms nicely under target 
space modular transformations, as expected for a string theory. It is not exactly modular 
covariant, but it turns out that F + is an anomalous modular form of weight 6g — 6. The 
deviations from modular covariance are of an interesting form, and they conspire to give 
an unexpectedly mild behavior at small area and/or weak gauge coupling. This may be 
related to a hidden symmetry in the string theory. The section begins with the derivation 
of a generating function for F+, followed by a general description of its modular properties. 
Then both are used to calculate F+ and F g for g = 1, ■ ■ • , 8, exactly in terms of modular 
forms. Finally, we examine the large and small area behavior of the free energy. 

QCD on the torus is simpler in many ways than on the sphere or on higher genus surfaces. 
The dimension factor drops out of the heat kernel expression for the partition function 



Z+ = Y, e-wE k KM-a)] e -»i ( 3 _ 22 ) 

7ii>7i 2 >...>n iv /2>0 

In the Gross- Taylor description this means that no omega points or omega-inverse points 
are required. The chiral partition function just counts simply branched maps. Despite the 
absence of omega points, the free energy is non-trivial since the map v : H g — > Sc =1 can 
wrap arbitrarily many times. The wrapping number n drops out of the Riemann-Hurwitz 
formula ( |2.20| ) for G — 1. This allows non-trivial modular transformations of the Kahler 
modulus A. In addition to formulations of QCD2 using gauge theory and string theory 
techniques that work for any Riemann surface, QCD on the torus has been reformulated as 
a two dimensional free fermion theory and a Jevicki-Sakita boson theory |9|, [10] . The relative 



simplicity of QCD on the torus along with the alternative formulations allows us to say a 
great deal about the string theory. 

The partition function ( |3.22| ) is similar to a theta function, as we saw above in (|2T8"|). It 
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would be nice to extract the large N expansion directly from the theta functions, but that 
has proved to be difficult. We will now develop a way to sum the series fl3.22| ) that is more 
amenable to a large N expansion. Since the sum is restricted to chiral representations, the 



factor of 1/N appearing in the first exponential in (|3.22|) is the string coupling. Z + may be 
expressed in terms of a generating function: 



Z+ 



J2 e ^Ej(A)(-\+ 1 - a )l e -^ Pknk 

7ii>7i 2 >...>n iv /2>0 



j= A 



V e ~^^( k < +k29a * +2kda * 527=k+i d "i) e -12« k m h 



m,->0, Vj 



(3.23) 



0.3=3 A 



where m k = n k - n k+ i (i.e. n k = £^ fe mz), j3 k = a t - a t -i and dp k = E£U#a r All of the 
sums over the rrij begin at zero, so the generating function is just a product of geometric 
sums. Summing these series (and allowing for a zero point energy eo) we find 



e e A/24 exp 



e^ A ' 2A exp 



XA 

"2N 



£ kdl + k 2 d ak + 2kd ak Y, d ai 



k=X 



l=k+l 
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XA °° ( 1 °° 

^T7 E k2d o k + 9 A 9a k - - E \ k - 1 \ d oA 

ZJV fc=l \ z 1=1 



Oj=jA 
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77(a) 



where the generalized eta function is given by 



3 A/24 



V(<*) y _j 



Sr 



+ 



e "^ 



Oj=jA 

(3.24) 



(3.25) 



The omitted terms correct for the fact that we have allowed an infinite number of rows instead 
of stopping j at N/2, the maximum for a chiral tableau. These terms are proportional to 
some power of e~ N , so they are exponentially small non-perturbative corrections. We can 
ignore them. When ( |3.25|) is evaluated at <x, = jA in the large N limit, the generalized eta 
function becomes the usual eta function 



V\ aj=j A — V 



q i/*Y[(i-<r) 



(3.26) 



n=l 



with q = e~ A (so that the Kahler modulus of the torus is r = — A/ (2m)). This is the genus 
g = 1 chiral partition function for SU(N) and U(N) Yang-Mills on the torus found in M. 

The free energy is the "connected" part of (|3.24|) . It may be computed using the 
Campbell-Baker-Hausdorff formula. This is particularly simple because log 77(a) = — ^ + 
yj,log(l — e aj ), and different derivatives cannot act on the same term in the sum. 
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It is key for a closed string formulation of QCD 2 that the odd powers of 1/iV vanish in 
the partition function. This leads to some non-trivial identities [fH|]. For example, the term 
in ( |3.24j ) proportional to 1/N (corresponding to genus g = 3/2) is 

^J k(k-l)q k 2kq 2k , ^ (k + l)q k + kl + l \_ n 

h\ W W* fe(l-? fc )(l-9')J ( } 

The argument given above for the absence of these terms, based on Gross's original discus- 
sion, is essentially group theoretic. It would be interesting to understand the identity ( |3.27| ) 
and its generalizations from the point of view of number theory. 

The free energy is almost, but not quite, invariant under modular transformations. Al- 
ready we have seen that F^~ = —log?]. This is not modular invariant, contrary to what 
one would expect for a string theory. Under the transformation of the Kahler modulus 
A — > 4tt 2 /A, the free energy has a modular anomaly: logr/(47r 2 /y4) = logrj(A) + 1 log(A/27r). 
It becomes modular invariant if we make an ansatz that the theory has a holomorphic 
anomalyQ The anomaly makes a contribution to F + that depends on A as well as A. Then 



Ft = -\ogN\{A + A)r 1 {A)r 1 {A)\. (3.28) 

This is modular invariant, and the original free energy is recovered by sending A — > oo, up to 
an infinite shift in the zero point energy. Modular invariance is an extremely useful property. 
It is a kind of strong/ weak coupling duality in the gauge coupling A, and it determines the 
general structure of the free energy in terms of modular forms. 

The higher genus free energies are also expected to be modular covariant after the non- 
holomorphic completion. The string coupling transforms non-trivially (with weight —1), so 
the higher genus free energies should transform with a definite modular weight. Douglas || 
was able to show that this is the case using his formulation of QCD 2 on the torus in terms 
of a Jevicki-Sakita two dimensional boson system: 

Z+ = f z>Lp e~ f dip ^ +9stA{dip)3 . (3.29) 

The propagator is {d<p(z) d<p(0)) = -d 2 logi^i) - ^ = p(z\A) + n 2 E 2 (A)/3 - |^, 
where $1 is the Jacobi theta function and the Weierstrass p function is given by jp(z|A) = 
z~ 2 + J2k 2(2& — 1) C(2k) E 2 k{A) z 2k ~ 2 . The functions E 2 k are the Eisenstein series-weight 2k 
modular forms. They are described below. The point is that the propagator has modular 



2 The holomorphic anomaly arises in topological string theory, coming from the contribution of BRST 
exact operators (with anti-holomorphic couplings) at the boundary of moduli space. |5| 
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weight two, where a weight k modular form transforms as Mfc(— 1/r) = r k Mk(T). This 
determines the weight of the genus g free energy 



1 fXA^ 29 ' 2 /r /■ „ i 2 5-2 
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(2(7-2)! V2iV; 



: (ty)< 



] . (3.30) 

/ connected 



It has 3(? — 3 propagators, giving a total weight 6g — 6. 

It would be fantastic if we could do the integrals to determine which weight 6g — Q 
modular form F + is. Unfortunately, these integrals are very difficult. Douglas has done 



the F 2 integral [15], and we have done the _F 3 integrals. The higher genus integrals are 
extremely difficult, since the convolution of offset p functions is not elliptic. Also, as the 
genus g increases, there are a growing number of diagrams and integrals. In fact, the number 
of diagrams at genus g is proportional to g\, as easily seen in the zero dimensional (p 3 integral. 
In addition, ever higher weight modular forms entering from the p function must be reduced. 
It does not seem promising to calculate the higher genus free energies this way. Fortunately, 
the fact that they are modular forms determines them up to a few coefficients. These 
coefficients may be found by taking parametric derivatives of the generalized rj function 
(pi). 



// 



A set of 3g 2 /4 ((3g 2 + 2g — 5)/4 for odd g) modular forms comprises a basis at weight 
6(7—6. The basis is made of products of E2, E4 and E§ or E2, E 2 ' and E 2 , which have weight 
2, 4 and 6, respectively. Most texts on modular forms write the basis in terms of E± and Eq, 
but for our purposes the second basis is preferable. The two bases are interchangeable, since 
E 4 = (E 2 f + 12E 2 ' and E 6 = E 2 3 + 18E 2 E 2 ' + 36E 2 ". Also, note that E 2 m = \[E 2 f - E 2 E 2 ' 
and Ek+2 = E^E-i + ^uE' k for k = 4, 6, • • •. The prime means differentiation with respect to 
A, so /' = d^f = j^d T f, and for finite A it becomes the covariant derivative 

D A = d A + ^rK:. (3.31) 

A + A 

acting on a weight k modular form. Regardless of which basis is used, the genus g free energy 
is determined up to roughly 3g 2 /4 coefficients 

f\A\ 2g ~ 2 3g - 3 
*?=(£) £ E c u E>{Etf(E») m . (3.32) 

\ ZiV / fc=0 2l+3m=3g-3-k 

The coefficients cu are rational numbers to be determined. 
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The k th Eisenstein series, Eu, is 



rjk ~ o/-fu\ ^ 



2A; ~ n k - 1 q n 2k ~ 

= 1 -5S^ :r = 1 "SS^ l(n)fl 

for fe G 2Z + (cf. fll6|). The number theoretic function akin) is the sum of the k th power 
of the divisors of n, <Jk{n) = J2d\nd h i an d -Bfc is the k th Bernoulli number. Every E\. is a 
modular form of weight k, except E 2 . There is no modular form of weight two, although E 2 
comes close. The Eisenstein series transform as 

/ 12 \ 

E 2 (-l/r) = t 2 (# 2 (r) + ^— j E k (-l/r)=r k E k (r) k = 4, 6, • ■ • (3.34) 

The fact that the higher Eisenstein series are modular forms is easily seen from the definition 

( [3.331 ) since the sums converge absolutely (and uniformly) on the upper half plane. Since E 2 

has a modular anomaly, it must be covariantized by adding a term depending on A 

oo 19 
E 2 = 24 d A F+ = 1 - 24 V <n(n) g n - - (3.35) 

n=l A + A 

so E 2 = 1 - 24e~ A - 72e- 2A . 

The G = 1 free energy for chiral U(N) (restricted to SU(N) tableaux) has been calculated 
exactly up to worldsheet genus 8. It would be possible to continue up to genus 11, but beyond 
that the computation takes too long even on fast computers. A glance at the genus 8 result 
( [A. 56 ) reveals that these expressions would fill pages. The first few free energies are given 



by 



F+ = ^A-log V (3.36) 

F+ = \ 2 f >- \hEl - 3E 2 E 4 - 2E 6 \ 

2!2434 V 2 ' (3.37) 



-XA 

2N 



2! 90 



{E 2 E 2 + E 2 } 



and 



(^Y 

Fs+ = ^^{- QE 2+ 1 5EiE 4 + 4E 2 3 E 6 -l2E 2 2 E 2 -l2E 2 E 4 E 6 + 7E? + 4E£} 



■i 



^{7[E 2 f + 3[E^f} 

(3.38) 
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Fi~ and F 2 + have been calculated previously by Gross @ and by Douglas [[H|, respectively. 
The higher free energies were not known. The remaining expressions F^, • • • , F 8 + are in- 
creasingly lengthy, so they are listed in the appendix. 

Using ( |2.15|) and ( p.!7| ) these results for the chiral U(N) free energies may be converted 
into the more complicated chiral SU(N) and full SU(N) free energies. The additional terms 
in the full SU(N) free energy coming from turning off the U(l) coupling and combining 
the chiral sectors are still modular forms (up to the E 2 modular anomaly), but they have 
lower weights. The corrections for the chiral SU(N) free energy do not have a definite weight 
because F +/ = (E 2 — l)/24 enters fl2.15|) . The same is true of the non-chiral U(N) free energy. 
The general form of the full SU(N) free energy (the analog of ( |3.32| ) for chiral U(N)) is 

g— 1 3g— 3— / 

F 9 = g s t 9 ~ 2 £(AA) 2 *- 2 -' J2 E 2 k M 6g _ 6 _ 2 ,_ 2M (3.39) 

1=0 k=0 

where Mk,i is a true modular form of weight fcora weight k combination of E 2 and E 2 in 
the other basis. For example, 



-\a \ 2 ( XA 



ijg- {7 [Etf + 3 |B»| 2 } - ^"I^U EiE, 



2^2 



F\ through F 4 are given in the appendix. The higher SU(N) free energies are omitted for 
brevity. 

A surprise about the free energy is that it has an unexpectedly mild singularity as XA — ► 0. 
To see this it will be useful to know the large and small area behavior of the Eisenstein series. 
The series are normalized so that 

E k -> 1 + 0{e~ A ) as A -»• +oo (3.41) 

for k — 2, 4, ■ • -. The small area limit is related to this by a modular transformation 

E k -> r- k + 0(e- 2ni/T ) as r -> i0 + (k = 4, 6, • • •) 



-► (-i/2vrz)- fc + C( e - 4 " 2/A ) as i^0 H 

1 1 9 

E 2 -► - - -— + 0(e~ 2m/T ) as r -> *0 + 

- =^ + S + 0(e" 4 ^) as i - 0+ 
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(3.42) 



It is amusing to note that these formulas also result from using the Euler-Maclaurin formula 
on the g-expansion ( |3.33| ). The usual derivation of the modular transformation laws is quite 



different. 

Consider the expression for the partition function fl3.24|) . In the large area limit _F + goes 



like e~ 2A for g > 2. There is no constant term, which would come from non-covering maps. 
Similarly there is no single cover term proportional to e~ A , since higher genus surfaces must 
be at least double covers of the torus. The absence of the constant term shows up in the 
expressions for F+ since each term contains a factor of E^ or E![ which begin with e~ A . It is 
part of the reason that the expressions are simpler in this basis. The small A limit is much 
more interesting. A weight 6g — 6 modular form would go like (1/A) 6g ~ 6 as A — > + . A 
much softer singularity is observed. Up to terms of order 0{e~ 1 ' A ) the following expressions 
hold: 

^ = f_^V(X-li + iil (3.43) 

2 V 2N J 13A 3 3A 4 45A 5 J V ; 

_/ XA\ 4 (-8 jW^IOOtt 4 224tt 6 1 
3 '" \~2NJ \I^ + ^47 9F + 81 A 9 J ^' ' 

F+ _ ( \A\ 6 (2272 2272 vr 2 8096 tt 4 41504 tt 6 48256 vr 8 1 
4 \~2NJ \VF~ 3 A 10 + 9 A 11 81 A 12 + 405 1 13 J 



„, / XA\ ( -13504 54016 vr 2 834304 tt 4 7010816 tt 6 
Ft = \-—\ { — + — ~^T— + 



2NJ I A 12 A 13 9A U 81 A 15 

17887904 vr 8 11958784 tt 10 1 
+ 



(3.46) 
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-+ 



405 A 16 1215 A 



\A\ 10 { 15465472 15465472 tt 2 105156608 tt 4 418657280 tt 6 



+ 



(3.47) 



V 2N J I 15 A 15 3A ie 9 A 17 27 A 18 

572409344 tt 8 2467804672 vr 10 33778284544 tt 12 1 
45 A 19 405 A 20 25515 A 21 J 

The expressions for Fj~ and _F 8 + are easily calculated as well from (|A.55|) and ( |A.56|) . 
In general the small area behavior of the free energy is 

/ \A\ 29 ~ 2 4g ~ 3 r, 7r2(fc-3(/+3) 

*-(-$ E,^^ (,48) 

where each Ck t9 is a rational number. This form may be proven using the Campbell-Baker- 
Hausdorff formula to get F+ from Z + ( |3.24 ), and then using the Euler-Maclaurin formula 
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to replace the sums with integrals. This small area behavior is interesting for a number 
of reasons. It might be the result of a symmetry in QCD 2 , or one in the underlying string 
theory. If it is, the symmetry would be novel. The small area limit is also interesting because 
of its implications for an equation relating the free energy at a given genus to that at lower 
genera. If the free energy satisfies a string master equation like the holomorphic anomaly 
equation, it must satisfy it as A — > 0. It is very easy to check that no simple equation will 
work. Of course, the equation could be more complicated. For example, if the full SU(N) 
partition function satisfies the simple holomorphic anomaly equation, then the chiral U(N) 
partition function will satisfy a non-polynomial (in 1/N) differential equation resulting from 



4 Conclusions 

In the recent program to extract a string theory from Migdal's explicit solution of QCD2, we 
have taken a step backward from the starting point, in a sense. We have computed even more 
explicit expressions for the QCD 2 string free energy up to genus 8. These calculations relied 
on a strong/weak gauge coupling duality that is exact at each order of string perturbation 
theory (but is violated by the non-perturbative corrections). The modular structure of the 
free energy is familiar from topological string theory, but there does not seem to be a simple 
holomorphic anomaly equation for F + . 

It might be expected that there would exist a handle generating operator since F + receives 
no contribution coming from the boundary of moduli space (collapsed handles or tubes). In 
some sense the differential operator generating the partition function from the generalized 
eta function ( |3.24| ) plays this role. It is not as simple as one would like, since it couples to 
each row number separately. This is equivalent to having it couple to the infinitely many 
deformations of QCD 2 (the higher Casimir perturbations), rather than coupling to the Kahler 
modulus A alone. So it is an important open question in the worldsheet theory to understand 
how the free energy at a given genus is related to that at lower genera. 

Another interesting question we have raised is the cause of the softening of the XA — » 
singularity. We proved the property by looking directly at the small area limit of the heat 
kernel expression for the partition function. It would be very interesting to have a worldsheet 
explanation for this effect (although it might just be accidental). 

In any case, the exact expressions for the free energy offer many possibilities for further 
investigations. The goal of the recent work on two dimensional Yang-Mills theory is to make 
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progress toward understanding four dimensional QCD, or at least to learn more about string 
theories without spacetime gravity. We have exhibited properties of the two dimensional 
free energy that could have a bearing on either of these two interesting goals. 
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A Appendix: The Free Energy up to Genus 8 

The U(N) free energy on the torus restricted to one chiral sector is 

F+ = -||i-log77 (A.49) 

F+ = -±*?y r {E 2 El + E»} (A.50) 

-XA^ 



F 3 + = ^-L {7 [E 2 f + 3 [£ 2 "] 2 } (A.51) 



-XA 



Ft = ±^ {27E 2 [E>] A - 36 [E 2 ] 2 [E>] 2 E>>- 746 [E'fE'i 
+ 12 [E 2 f [E 2 "] 2 + 2A6E 2 E 2 ' [E 2 "\ 2 - 106 [E 2 n f] 

(^XA\* 

F+ = ^H^ {162 [£ 2 ] 4 [£ 2 f + 5940 [E 2 ] 2 [E 2 f + 85299 [E 2 f + 19821 [E 2 "} 4 



- 216 [E 2 f [E 2 X E 2 ' - 8352 [E 2 Y [E 2 '} 6 E' 2 ' - 165388£ 2 [E 2 f E' 2 ' 
+ 72 \E 2 f [E 2 n \ 2 + 2928 [E 2 ] 4 E' 2 [E 2 n ] 2 + 97980 [E 2 f [E 2 'f [E 2 "] 2 
+ 291334 [E 2 f [E 2 "} 2 - 16896 [E 2 f [E 2 'f - 116472 [E 2 ] E' 2 [E 2 f] 



(A.53) 
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-aW° 



1 v J {216 [E 2 ] 7 [E 2 f + 12420 [E 2 f [E 2 f + 935442 [E 2 f [Ef 



in| 07 ^ iJJ -" L-^^J L"^2J ' — — " L-^J L^2J ' """■■■■■■■" L-^^J L W 2J 

+ 9054978^2 [£ 2 'f - 288 [E 2 f [E 2 f E 2 " - 17136 [E 2 f [E 2 f E 2 " 

- 1798164 [E 2 f [E 2 f E' 2 - 21957340 [E 2 f [E 2 f E 2 " - 49812944 [E 2 'f E 2 
+ 96 [E 2 f {E 2 "] 2 + 5904 [E 2 ] 7 E' 2 [E 2 n \ 2 + 1129560 [E 2 f [E 2 f [E 2 'f 
+ 16998104 [E 2 f [E 2 f [Elf] 2 + 85070724 [E 2 ] [E 2 f [E 2 f - 230768 [E 2 f [E 2 n f 

- 4242768 [E 2 f E 2 [E 2 "f - 43083696 [E 2 ] 2 [E 2 f [E 2 "] 3 - 55574424 [E 2 f [E 2 f 

+ 5750160 [E 2 f \E 2 n f + 22892460^2^' {E 2 "} 4 - 2132916 [E 2 "} 5 } 

(A.54) 



-A^ 12 



jtf 



F+ = \ 2 ^ y (2592 [£ 2 ] 10 [E 2 f + 204768 [£ 2 ] 8 [E 2 f + 89976744 [£ 2 ] 6 [£. 
12! 81 *■ 

+ 1909213524 [E 2 f [E 2 ] 7 - 3456 [E 2 ] u [E 2 f E' 2 - 279936 [E 2 f [E 2 'f E' 2 ' 

- 177698880 [E 2 ] 7 [E 2 f E' 2 ' - 4174329312 [E 2 f [E 2 f E' 2 ' 

- 62394892148 [E 2 ] 3 [E 2 f E' 2 ' + 1152 [E 2 ] 12 [E 2 "] 2 

+ 95616 [E 2 ] w E 2 [E 2 f + 116614944 [E 2 ] 8 [E 2 ] 2 [E 2 "] 2 

+ 3013923216 [E 2 f [E 2 f [E 2 n ] 2 + 63439748328 [E 2 ] 4 [E 2 f [E 2 n ] 2 

- 25422720 [E 2 f [E 2 n f - 719490240 [E 2 ] 7 E' 2 [E 2 "] 3 

- 26470959696 [E 2 f [E 2 f [E 2 n \ 3 - 216898783824 [E 2 f [E 2 'f [E. 
+ 3707598864 [E 2 f [E 2 f + 44962956000 [E 2 ] A E 2 [E 2 n f 

- 27342815040 [E 2 f [E 2 f + 21435613473 [E 2 ] 2 [E 2 f 
+ 354642271752 [E 2 ] 2 [E 2 f [E 2 n ] 2 + 246384220752 [E 2 ] 2 [E 2 f [E 2 n f 

- 221060144076^2 [E^] 7 E' 2 - 540228415584£ 2 [E 2 f [E. 

- 75540784608^2^2 [E 2 f + 37289952912 [E 2 f 
+ 343299239380 [E 2 f [E 2 "] 2 + 186490620756 [E 2 f [E. 

+ 4465217052 [E 2 'f} 



■J"] 3 



^ 3 



*,T 



(A.55) 



-XA^ 14 



F+ = V ^ {31104 [E 2 ] 13 [£ 2 '] 4 + 3146688 [£ 2 ] n [£ 2 f + 9381744000 [£ 2 ] 9 [£ 2 f 

+ 284491245600 [E 2 ] 7 [E 2 ] 7 - 41472 [E 2 } u [E 2 f E' 2 ' - 4278528 [E 2 ] 12 [E 2 f E' 2 

- 27941270346864 [E 2 f [E 2 f E 2 - 606326798592 [E 2 ] 8 [E 2 f E' 2 

- 18712973376 [E 2 ] 10 [E^f E' 2 ' + 13824 [E 2 ] 15 [£ 2 'f + 1453824 [E 2 ] 13 E' 2 [E 2 "] 2 
+ 12436238592 [E 2 ] u [E 2 f [E 2 f + 428976125568 [E 2 f [E 2 f [E 2 "f 

+ 27464636483952 [E 2 ] 7 [E 2 f [E 2 "} 2 - 2753687808 [E 2 ] 12 [E%f 

- 100795663872 [E 2 ] 10 E 2 [E 2 n f - 11670697945728 [E 2 ] 8 [E 2 f [E 2 n f 

- 191631253050048 \E 2 f [E 2 f [Egf + 1796319251712 [E 2 ] 9 [E 2 f 
+ 35809340742912 [E 2 ] 7 E' 2 [E 2 "] 4 - 65393019413760 \E 2 f [£ 2 "] 5 

+ 10411506055692 \E 2 f [E 2 f + 374290947363144 \E 2 f \E 2 f [E 2 "] 2 

+ 622325551111920 [E 2 f [E 2 ] 2 [E 2 f - 314884801430312 [E 2 ] 4 [E 2 ] 7 E 2 ' 

- 1915319874231576 [E 2 ] 4 [£ 2 '] 4 \E 2 n f - 599386224758688 [£ 2 ] 4 E' 2 \E 2 n f 
+ 95567800131858 \E 2 f [E 2 ] 9 + 2553509520593988 [E 2 f [E^f [E 2 f 

+ 3327996287536944 [E 2 f [E^f [£ 2 "] 4 + 207585191777040 [E 2 f [£ 2 "f 

- 1483200593325666 [E 2 f [E^f E' 2 ' - 6326492920454424 [E 2 f [£ 2 '] 5 [£ 2 "f 

- 2174062501379952 [E 2 f \E 2 f [E 2 f + 287886256181076£ 2 [E 2 f° 

+ 4807697809119300E 2 [E 2 f ] 7 [E 2 "} 2 + 5189221094011812£ 2 [£ 2 '] 4 [£ 2 "] 4 
+ 436733614643256£ 2 £ 2 [E 2 n f - 1172778987525768 [E 2 f E 2 " 

- 3523420607226032 [E 2 f [E 2 f - 1116320466888120 [E 2 f [E 2 f 

- 17870538853512 [£ 2 "] 7 } 

(A.56) 
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The full SU(N) free energy on the torus is 



F 1 = ^(^)-21ogr ] (A.57) 



{£ 2 £ 2 ' + E^} + ^_Ze 2 ' (A.58) 



2! 45 L z /J 2! 6 



4! 27 1 J 2 4!3 (A.59) 



H^{4£ 2 £ 2 ''-5[£ 2 '] 2 } 



+ 
+ 
+ 



F4 = ^27 ^ [E ^ ~ 36 ^ [ ^ ]2 E * ~ 746 ^ ^ 2 " 

+ 12 [E 2 f [E 2 n f + 2AQE 2 E 2 > [E 2 n f - 106 [£ 2 "] 3 } 

(=mY(m.) 

K2N \ K2m) f 150 [£'] 4 - 160£ 2 [£'l 2 £'' + 40 \E 2 f [E'f + 180 E' [E'f] 
6! 3 L " J 

V2jvy , V2JV2y (30£ 2 [£ 2 'f - 20 [£ 2 ] 2 E'E 2 " ~ 120 [£ 2 '] 2 £ 2 " + 60E 2 \E 2 "] 2 } 
6! *• . - j 

^-f {65 [E 2 f - 360 [£ 2 ] 2 [E 2 f + 240 [£ 2 ] 3 £ 2 " + 420£ 2 £ 2 '£ 2 " - 480 [E 2 n f } 
6! 24 l > 

(A.60) 
The free energy for genus 5, 6, 7 and 8 may be computed from ( |A.53|) to ( |A.56j ), but they 
are omitted to save space. 
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